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We investigate the electrical transport properties of two hydrogen tautomer configurations of 
phthalocyanine (H2PC) connected to cumulene and gold leads. Hydrogen tautomerization affects 
the electronic state of H2PC by switching the character of molecular orbitals with the same symmetry 
close to the Fermi level. The near degeneracy between the HOMO and HOMO-1 leads to pronounced 
interference effects, causing a large change in current for the two tautomer configuratons, especially 
in the low-bias regime. Two types of planar junctions are considered: cumulene-H2Pc-cumulene and 
gold-H2Pc-gold. Both demonstrate prominent difference in molecular conductance between ON and 
OFF states. In addition, junctions with gold leads show pronounced negative differential resistance 
(NDR) at high bias voltage, as well as weak NDR at intermediate bias. 
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I. INTRODUCTION 

The ability to control current in molecular junctions 
is a prerequisite for functional devices in molecular 
electronics. Of particular interest for nanoelectronics 
are two-terminal switches characterized by two distinct 
current- volt age characteristics (JVC) representing ON- 
and OFF-states, switchable by exceeding threshold volt- 
ages in hysteretic junctions. With such two-terminal 
components one can in principle build nanoscale mem- 
ory and resistor-diode logic for scalable architect ureJ^^. 

A major challenge is to develop two- terminal unimolec- 
ular switches (UMS) where the switching is truly in- 
tramolecular, and not due to cooperative or interface 
effects in multi-molecule junctions. One UMS example 
is provided by a recent experiment by Lortscher et al.^ 
using a mechanically controlled break junction (MCBJ) 
to demonstrate switching and hysteretic current-voltage 
characteristics (IVC) of a single Tour- wire molecule^ ^. 

An UMS should be electrically coupled and immobi- 
lized in the circuit, yet have enough freedom to perform 
switching. Up to now, several different kinds of molecu- 
lar switches have been proposed, most of them involving 
a change of molecular conformatioiP^, or redox stat#^. 
An interesting case related to the present discussion is 
an SnPc switclP where the current is determined by the 
vertical position of the tin atom, controlled by a scanning 
tunneling microscope (STM). 

Recently, Liljeroth et al^^ proposed a new type of 
molecular switch where the switching is mediated by a 
tautomerization reaction of metal- free naphthalocyanine. 
Using a low-temperature STM, a voltage pulse at the 
STM tip can induce a change in the orientation of the 
hydrogen atom-pair at the center of naphthalocyanine, 
leading to switching between low and high conductance. 
The fact that naphthalocyanine does not change posi- 



tion or conformation of the molecular framework could 
be very beneficial for controlling molecular electronic cir- 
cuits. If the molecule can be integrated in the circuit, it 
will not move or change its outer shape upon switching. 

There is now a variety of techniques that can be used 
for contacting and investigating single molecules (see e.g. 
Refs.im^). However, when it comes to interconnecting 
single molecules to functional electronic circuits, there 
are currently no working solutions. The circuits will nec- 
essarily have to b e base d on various forms of self-assembly 
of building block^^^, but the field is wide open, and ex- 
tensive fundamental research will be needed before useful 
functional systems can emerge. 

One possibility to wire circuits could be to first deposit 
the molecules and then grow the wires. This case has re- 
cently been investigated by Okawa et a/.^^^K^, with the 
pronounced goal to fabricate a single-molecule integrated 
circuit. Nano-clusters of phthalocyanine molecules are 
adsorbed on a molecular layer of diacetylene, and chain 
polymerization is initiated by applying a pulsed bias volt- 
age to the row of diacetylene molecules to which a ph- 
thalocyanine pentamer is adsorbed. In this way Okawa 
et al¥^ recently succeeded in connecting two polydi- 
acetylene chains to the same pentamer phthalocyanine 
molecule. 

In this paper, we focus on the planar electrical trans- 
port properties of a UMS based on hydrogen tautomer- 
ization in metal- free phthalocyanine (H2PC) connected 
to semi-infinite metallic wires of carbon (cumulene^P^ or 
gold. The transport properties are examined as a func- 
tion of the electronic structure of the tautomer state of 
the molecule, the type of leads, and the coupling be- 
tween electrodes and the molecule. Our results show that 
switching the orientation of the H-H atom pair in the 
H2PC cavity can significantly change the conductance of 
the molecule, effectively leading to ON- and OFF-states 
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for both cumulene (Fig. [T]) and gold (Fig. [2| leads. More- 
over, in the case of gold we find that negative differential 
resistance (NDR) appears at high bias voltages due to 
reduced transmission through one of the contacts. 

A few important clear limitations should be noted from 
the start: we treat neither (i) the influence of any sub- 
strate nor (ii) the switching mechanism and stability 
(usefulness) of the switch. Liljeroth et al}^ have demon- 
strated experimentally that the switch works in principle 
in vertical transport (STM) with the molecule lying flat 
on an insulating monolayer on top of a metallic substrate. 
In our work we instead probe theoretically the horizon- 
tal (planar) transport properties of a "freestanding" 2D 
system of a H2PC molecule plus leads. Connecting molec- 
ular devices in the 2D plane will be essential for building 
functional circuits - here we present some first results for 
planar transport of some simple model systems in order 
to assess their potential usefulness as switches or non- 
linear elements. 

The mechanism of bias-driven tautomerization h as re- 
cently been attributed to electron-phonon couplin^^^^, 
implying that the switching can most probably be trig- 
gered by either an STM-like current pulse into the molec- 
ular plane, or a current pulse along the H2PC plane. 
The results of Liljeroth et al}^ showed the STM-driven 
molecule to be bistable, however fluctuating (telegraph 
type of noise) between the two conductance levels with 
time (and especially under elevated bias). We would like 
to point out that the interaction of the H2PC with the un- 
derlying substrate stabilizes the OFF and ON states by 
increasing the transition barrier. In the related bistable 
case of Ag-Ceo-Ag^'', the hysteretic behaviour was more 
stable, except at high bias where fluctuating switching 
dominated. 

As stated at the beginning, a switch has two prop- 
erties: (i) an IVC with two branches, ON and OFF; 
and (ii) a mechanism for switching between those two 
branches. The OFF state can typically be generated by 
(a) interrupting a single conduction path by disrupting 
wave function overlap; or (b) by interference of two con- 
duction paths, placing an interference minimum at the 
Fermi level. 

In the present calculation, the difference between the 
ON and OFF states is produced by a two-channel inter- 
ference effect, involving nearly degenerate HOMO and 
HOMO-1 orbitals and controlled by the H-H tautomer 
switch. In a molecular context, similar interferen ce ef- 
fects were discussed early on by Sautet and JoachinJ^^'^ 
for transmission through benzene described as two-path 
interferometer, and has recently been elaborated for ben- 
zene and other systems in several papers using state of 
the art methods^^ In all of these cases, the basic idea 
is to influence the interference by varying the phase shift 
along one of the paths, controlling the transmission at 
the Fermi level and implementing a transistor or switch. 

In the present paper we study the case where the two 
interfering states remain pinned at the edge of the bias 
window. Therefore the entire IVC is controlled by in- 



terference: the zero-bias destructive interference mini- 
mum in the OFF-state stays within the voltage window, 
following the Fermi level of one of the leads, while the 
transmission peaks stay outside. In the absence of the 
other transmission peaks (conducting states) in the volt- 
age window this makes it possible, in principle, to fol- 
low e.g. the OFF-state with increasing bias until the 
molecule switches to the ON-state, generating a hys- 
teretic IVCP^. Other methods for switching would be 
via current pulsed or photoexcitatioiP^^. 



II. COMPUTATIONAL METHOD 

To model the II2PC junctions, in principle the external 
leads are coupled directly to the outer rim of the H2PC 
molecule. However, the computational method requires 
that part of the leads are treated as an extension of the 
molecule, forming an extended molecule and scattering 
region coupled to semi-infinite leads. For coupling the 
semi-infinite cumulene leads to the H2PC molecule (Fig. 
[1]) we attach two cumulene-like spacers -CH=C=C=C via 
a TT-bond, expected to provide a stable and transparent 
coupling. For coupling semi-infinite gold wires to H2PC 
(Fig. [2|, we model this situation by using short atomic 
Au chains, Aun (n=l and 2), connected to Au(OOl) leads, 
as well as one case of semi-infinite Au chains. 

H2PC is connected to two different kinds of electrodes, 
carbon and gold chains. Such two-probe systems are di- 
vided into three regions: left and right electrodes and 
scattering region. The scattering region contains a por- 
tion of semi-infinite electrodes to account for electronic 
and geometric relaxation at the metal-molecule interface. 
The geometrical setup procedure of two electrode mod- 
els and computational details can be briefly described as 
follows: 

(i) cumulene-H2Pc-cumulene 

First, the free H2PC molecule was optimized. Two iso- 
mers (a) H-H pair horizontal (H2Pch) (b) H-H pair ver- 
tical (H2Pcv) were considered. To setup the cumulene 
chain electrode system, one hydrogen atom per contact 
site of optimized H2Pch and H2PCV was removed and sub- 
stituted with a (-CH=C=C=C) unit on each side, char- 
acterized as a scattering region of the extended molecule. 
The substituent unit was linked to H2PC at the peripheral 
group. All configurations, all atoms of scattering region 
were relaxed. The scattering region was sandwiched be- 
tween semi-infinite carbon linear chains, cumulenes, con- 
sisting of carbon double bonds. Optimum molecule-lead 
distances were obtained via energy minimization. 

(ii) AU-H2PC-AU 

Possible adsorption configurations of the scattering re- 
gion were selected to closely resemble the experimental 
arrangement of Nazin et al^, where both short atomic 
gold chains and the molecule were assembled on a metal- 
lic NiAl(llO) substrate. In our model, the gold atoms 
were placed at 77^ sites of opposite benzene rings of H2PC. 
We considered different models for the gold leads: (a) 
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Au(OOl) semi-infinite needle leads connected via short 
gold chain, labelled as H2Pch/H2Pcv@Aun-Au(001), n=l 
and 2, (b) semi-infinite gold chain lead coupled to the 
H2PC molecule, labeled as H2Pch/H2Pcv@Au-chain. For 
all systems, the Au-Au positions and bond length of the 
gold atomic chain in both scattering region and the leads 
were fixed at 2.88 A^^. The short Au chains were con- 
nected to the Au(OOl) electrodes at the hollow site. 

The geometry optimization of all systems above has 
been implemented in the localized-atomic-orbital code 
SIESTA package^^ The GGA-PBE^^ and LDA-PW9240 
approximations were used for the exchange correla- 
tion functional of cumulenes and gold electrodes, re- 
spectively, and real-space integrations was performed 
using 200 Ry cutoff. The atomic core electrons 
were modeled with Troullier-Martins norm-conserving 
pseudopotentiaPH, and the valence states were 2s2p for C 
and N, Is for H, and 5d6s for Au. In both models, the su- 
per cells were created with large enough vacuum regions 
to prevent interaction due to mirror image. We used a 
double zeta basis with polarization orbitals (DZP) for 
H, C, N atoms and a single zeta basis with polarization 
(SZP) for Au atoms. On each atom, which was allowed to 
move, was relaxed until the force becomes less than 0.01 
eV/ A. Convergence of total energy with respect to k- 
points sampling was tested up to 3x3x1 Monkhorst-Pack 
k-point mesh, yielding 6 k-points in the BZ. 

To investigate electrical transport properties of two- 
probe systems, we performed non-equilibrium Green's 
function (NEGF) technique based on density functional 
theory-based theory (DFT) as implemented in the Tran- 
SIESTA package^^. The exchange-correlation potential 
was approximated by GGA-PBE and LDA-PW92 for cu- 
mulenes and gold electrode, respectively. DZP was used 
for H, C, N atoms while SZP was used for Au atom. 
When a bias voltage was applied to the junction, the 
electrical current was calculated by integrating the trans- 
mission function over the bias window. In this work, 
the current- volt age characteristics have been calculated 
in the bias region from -2 to 2 V, where the potential bias 
was incremented in steps of 0.2 V. 

In the present calculations, the HOMO-LUMO gap 
is about 1.2 eV (Figs. [3] and [4|, and the Fermi level 
lies in the vicinity of the HOMO, determined by self- 
consistent charge transfer. In the experimental setup of 
Liljeroth et al}^ the HOMO-LUMO gap is about 2.5 eV, 
while the gas phase gap is about 4 eV. It is well-known 
and understood that local approximations of the density 
functional underestimates the HOMO-LUMO gap, and 
consequently often overestimates the conductance. Nev- 
ertheless, NEGF-DFT calculations with standard local 
exchange-correlation potentials are widely used and pro- 
vide benchmark results, and we believe that the present 
calculations provide good estimates of current-voltage 
characteristics (IVC), especially the differences induced 
by switching between the tautomer states. A recent dis- 
cussion of HOMO-LUMO gaps in CuPc can be found 
in Ref. |43l and the results of more advanced (and time 



consuming) calculation using the so-called GW approxi- 
mation for benzene on graphite (0001) is given in Ref. [44l 
and for Au-S-benzene-S-Au junctions in Ref. l45l and l46l 



III. RESULTS AND DISCUSSION 
A. Cumulene-H2Pc-Cumulene 

Fig. [1] exhibits the I-V characteristics (IVC) and cor- 
responding differential conductance for the two cases of 
hydrogen atoms orientation, H2Pch and H2PCV, in the 
bias region from -2 to 2 V. The mirror symmetric struc- 
ture of the H2Pch and H2PCV molecules results in sym- 
metric I-V curves. H2PCV performs as the ON-state, and 
H2Pch as the OFF-state. 

We first discuss the T{E^ V = 0) zero-bias part of the 
transmission spectrum of H2Pch and H2PCV in Fig. [3] 
and Fig. |4]in some detail. The molecular projected self- 
consistent Hamiltonian (MPSH) eigenvalues (a-b in the 
OFF and c-d in the ON states, marked as lines in Fig. [3| 
correspond well to the transmission resonance peaks. 

The orientations of the molecular eigen-orbitals de- 
pend on the position of the pair of hydrogen atoms. As 
seen in Fig. |3j the rotation of the H-H pair switches the 
character of the HOMO and HOMO-1 eigenstates. The 
HOMO of the OFF-state (a) is weakly coupled to left 
and right leads and gives rise to a narrow transmission 
peak close to the Fermi level (^=0), while the HOMO-1 
(b) is strongly coupled to the leads, resulting in a broad 
transmission peak well below the Fermi level (Fig. ^p). 

Turning the H-H pair knob to the ON-state reverses 
the orbital character of the HOMO and HOMO-1 lev- 
els and changes the orbital character at the Fermi level. 
Now, the HOMO of the ON-state (d) is strongly coupled 
to the left and right leads, resulting in a broad transmis- 
sion peak spanning the Fermi level (Fig. |4|l). 

The bias dependence of the transmission spectra 
T{E^ V) illustrated in Fig. [ijexplains the basic reason for 
the large ON/OFF conductance ratio in the 0-1 V bias 
range: the OFF-state transmission resonance stays al- 
most entirely below the Fermi level of the lead for which 
chemical potential /i = EF-eV/2 goes down under ap- 
plied bias, while the ON-state resonance extends into the 
bias voltage window. Apart from this, the transmission 
spectra are rather similar. 

Since the transmission near the Fermi level proceeds 
via two channels (orbitals), one can expect interference 
effects to affect line shapes and widths of the spec- 
tral peaks. Fig. |4] indeed provides a nice illustration 
of Fano resonances and interference in resonance trans- 
mission involving a narrow and a broad (continuum) 
resonanc^^^MU^ showing the whole range of lineshapes 
from a pure window resonance (interference minimum) ( 
Fig. Qa), V = 0) to strongly asymmetric resonances with 
pronounced interference minima (Fig. [4]^a), V > 0.4) 
(also seen in the high- lying LUMO peaks). More im- 
portantly, however, the interference has the most pro- 
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Bias voltage (V) Bias voltage (V) 

FIG. 1: (Color online) The calculated I-V characteristics of cumulene(E'F+eV/2)-H2Pc-cumulene(E'F-eV/2) junctions for two 
cases of inner hydrogen H-H pair orientation in H2PC, obtained in the bias region from -2 to 2 V. Note that the IVC is ohmic - 
this is not a tunnel junction. The right panel shows the conductance of ON-state and OFF-state and the ON/OFF conductance 
ratio. For V > 0, the conductance is calculated as I / V. The zero-bias conductance is given by the ratio of the ON and OFF 
transmission coefficients, and takes the value of 431. 



nounced effect on the spectral density at the Fermi level 
in the OFF-configuration, the Fermi level falling in an 
interference minimum at zero bias (Fig. gb), V = 0). 

The H2PC molecule therefore functions as an interfer- 
ometer where, by rotating the H-H pair by ninety de- 
grees, we reverse the asymmetry of the line shape of 
the transmission through the HOMO/HOMO-1 orbitals, 
placing an interference minimum well below the Fermi 
level (ON state), or in the HOMO-LUMO gap, suppress- 
ing the transmission at the Fermi level (OFF state). 

It should be noted that Sautet and Joachim de- 
scribed electronic interference produced by a benzene 
molecule embedded in a polyacetylene chain^^ and ap- 
plied this to a molecular switch controlled by photoin- 
duced tautomerization^^. The results demonstrated the 
sensitivity of the zero-bias transmission to the relative 
position and interaction of the states, disturbed through 
the tautomerization. Those interferometer effects are 
similar to the present case. An important difference, 
however, is the behaviour of such an interferometer under 
applied finite bias voltage. 

From the bias voltage dependence of the transmission 
peaks in Fig. [4] we can draw a number of conclusions 
regarding charge transfer, level pinning, and bias voltage 
drop over the junction. Starting with the ON-state in 
Fig. [4^, we note that the entire transmission spectrum 
shifts almost uniformly to lower energy in the direction 
of the downgoing lead with chemical potential ja = E^- 



eV /2. At zero bias ( F=0), charge transfer between H2PC 
and both leads places the HOMO resonance at the com- 
mon Fermi level (E'f=0) of the junction. With increas- 
ing bias F>0, Fig. |4^ shows that the HOMO-resonance 
becomes roughly pinned to the lead with downgoing po- 
tential, reflecting bias-dependent electron transfer from 
H2PC to the lead (positively charging the molecule) to 
keep the molecular level of H2PC equal to the chemical 
potential of the lead, /i = EF-eV/2. 

The bias voltage drop in Fig. [6^ confirms this picture: 
the main bias voltage drop occurs in H2PC near the left 
contact (top panel), and most of the bias voltage increase 
is distributed over the bulk of the H2PC molecule, the 
levels therefore following the chemical potential of the 
right, downgoing lead. 

Examining the ON state on Fig. [4^, one notes that the 
broad HOMO, the HOMO-2 and LUMO move according 
to -0.4V, not -0.5V, while the narrow HOMO-1 
and the high-lying LUMOs move as -0.5V. This re- 
flects the localization of the corrsponding orbitals: e.g. 
the HOMO-l(Fig. (sjc)) is well localized m a region con- 
trolled by the right lead (Fig. |6|, while the HOMO (Fig. 
|3jd)) extends into a region with upgoing potential /i = 
E'F + eF/2 controlled by the left lead. 

Now considering the OFF-state in Fig. (4)3, the motion 
of spectral peaks is the same as for the ON-state in the 
0-lV bias voltage region. However, above IV the lev- 
els begin to float, becoming independent of bias voltage. 



5 




-2 -1 1 2 0.2 0.6 1.0 1.4 1.8 

Bias voltage (V) Bias voltage (V) 



FIG. 2: (Color online) I-V characteristics for three cases of Au-H2Pc-Au junctions for two different H-H pair orientation. 
Note that the JVC is ohmic - this is not a tunnel junction. The right panel shows the conductance of H2Pc7;@Au2-Au(100) 
(ON-state) and H2Pc/i@Au2-Au(100) (OFF-state), and ON-OFF conductance ratio. For V > 0, the conductance is calculated 
as I / V. The zero-bias conductance is given by the ratio of the ON and OFF transmission coefficients, and takes the value of 
86. Also note the appearance of negative differential resistance (NDR) around 1.4 V bias in all cases, as well as weak NDR 
around 0.7 V in Au(001)-Au2-H2Pc-Au2-Au(001). 
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FIG. 3: (Color online) The transmission spectra, the MPSH 
eigenvalues and the MPSH eigenstates around Fermi level 
of the 'OFF' and 'ON' configurations of cumulene-H2Pc- 
cumulene. 



This signals that there is no longer any net charge trans- 
fer from H2PC to the right lead. Indeed, as seen in Fig. 
[gJd, the electrostatic potential does not change much over 
the central- left part of the molecule in the 1.6 < V < 2.0 
interval. 
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FIG. 4: Bias dependent transmission of cumulene-H2Pc- 
cumulene as a function of energy E. (a) ON-state; (b) OFF- 
state. The Fermi energy is set to zero; the shaded area indi- 
cates the bias voltage window. Notable features are (i) the 
transmission lines moving down, pinned to the downgoing 
right lead; (ii) the decoupling from the leads in (b) above 
1.2V bias; (iii) the two interfering transmission channel and 
the Fano lineshapes. 
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A closer look at Fig. [6^,b reveals some interesting sim- 
ilarities and differences in the voltage drop and screening 
between the ON- and OFF-states. In the ON-state (Fig. 
[6|i), the bias voltage drop is quite stable, but it neverthe- 
less shifts from the center to the left edge of the 5-ring 
pyrrole group in the < V < 2 bias interval. In the 
OFF-state (Fig. [gJd), the bias voltage drop is first lo- 
cated around the left edge of the 5-ring pyrrole group in 
the <V < 1 bias interval, and then jumps towards the 
center of the molecule in the 1.2 < V < 1.6 bias interval 
(and in fact shifts back a little in the 1.6 < F < 2 range). 

In the OFF-case the conduction state is localized over 
the central part of the H2PC molecule (essentially Fig. 
[sj^top, left)) and the voltage drop moves from the left 
edge to the center of the charge distribution. It is then 
understandable that the spectral peaks first follow the 
right lead and then become bias-independent. 

In the ON-case, on the other hand, the conduction 
state is delocalized over the H2PC molecule and the leads 
(Fig. |3|top, right)), and the bias voltage drop basically 
stays put at the center of the left part of the charge distri- 
bution. It is then understandable that the spectral peaks 
follow the right lead over the entire bias region. 

The HOMO coupled to the electrodes in the ON state 
provides an open system with good conductance and 
metallic-like screening. It is charge-transfer connected 
to the right lead and (de) charges to stay close to the 
chemical potential of the right lead. This picture does 
not change in the investigated < V < 2 range. 

On the other hand, the HOMO in the OFF state is a 
more closed system with less transparent coupling to con- 
tacts. As a consequence, the charge transfer to the right 
lead will be more difficult to maintain with increasing 
bias voltage. This means that the molecule will charge 
up (accept electrons) and float up in energy. This will 
shift the bias voltage drop toward the right electrode, 
as already discussed in connection with Fig. (6)3. These 
conclusions are, of course, dependent on the nature of 
the electrodes used to wire the molecule. To be able to 
identify the effect of the molecular-electrode coupling, we 
therefore now also consider gold electrodes. 



B. Au-HsPc-Au 

Fig. [2] presents IVCs for the two H2Pch/H2Pcv tau- 
tomer systems connected to several types of gold leads: 
3-D semi-infinite needles connected via short 1-D atomic 
Au-chains (the chains arranged as in the experimental 
setup by Nazin et al.^^)^ as well as 1-D Au wires. In 
the case of gold leads, the ON and OFF states corre- 
spond to the H2Pch and H2PCV tautomers, respectively. 
In comparison with cumulene, this reversal of ON and 
OFF states, relative to the orientation of the H-H pair, 
is due to the symmetry of the coupling orbitals at the 
H2PC-AU contacts. 

The detailed structure of the IVCs in Fig. [2] can be un- 
derstood from the bias- voltage dependence of the trans- 



mission spectra of Au(100)-Au3-H2Pc-Au3-Au(100), as 
shown in Fig. [sja) and Fig. [sjb) respectively. We first 
discuss the T{E, V = 0) zero-bias part of the transmis- 
sion spectrum in Fig. [5] in some detail. The linewidths 
with Au leads (~ 0.2eV) are considerably narrower than 
with cumulene leads (~ 0.4eF), manifesting a weaker 
coupling to the leads. It should be understood that the 
molecule is only physisorbed at the electrodes, bridging 
the gold chains arranged at the surfacd^. The hydro- 
gen tautomerization switches the orbital character of the 
HOMO and HOMO-1 levels. The profiles of the conduc- 
tion peaks result from the interference between HOMO 
and HOMO-1, giving rise to one symmetric broad peak 
with a dip (ON-state, much like in the case of the cu- 
mulene electrodes), and two separate peaks with steep 
edges (OFF-state). The spectrum around the Fermi level 
(-0.5 < £^ < 0) therefore derives from HOMO and 
HOMO-1, with different detailed orbital character upon 
hydrogen tautomerization, turning the H-H "knob". 

The zero bias MPSH HOMO and HOMO-1 orbitals of 
two tautomer systems are shown in Fig. [5] Similar to 
the cumulene wire leads in Fig. |3j the inner-hydrogen 
tautomerization of H2PC switches the molecular orbital 
character of HOMO and HOMO-1. One molecular or- 
bital shows the interaction between the gold d^j^-orbital 
and the carbon p-orbital (weak overlap), while the other 
one shows the interaction between the gold s-orbital and 
the carbon p-orbital (strong overlap). Concerning the 
transmission at the Fermi energy, in the ON-state the 
electrons can propagate from the left to the right elec- 
trode. In the OFF-state, the electrons are almost totally 
reflected (at the contact to the right lead, as can be shown 
with transmission eigenchannel analysis), resulting in the 
conduction being almost blocked. 

Again, from the bias voltage dependence of the 
T{E^ V) transmission peaks in Fig. [5] we can draw a 
number of conclusions regarding charge transfer, level 
pinning, and bias voltage drop over the junction, simi- 
lar to the case of cumulene wire leads: For H2Pch (ON- 
spectrum), at low bias voltage in the range < V < 0.4, 
the current increases because the resonance peaks move 
distinctly slower than ja = Eip-eV/2^ the bias window 
therefore including an increasing part of the line wings. 
In the range 0.4 <V<1, the transmission peaks move as 
fast as the chemical potential of the right lead, resulting 
in the plateau in the IVCs. For V > 1, something inter- 
esting happens: the intense resonance peaks shift back 
into the bias window and start following the chemical 
potential /i = E-p-\-eV/2 of the left lead, resulting in a 
large increase of the current in the range 1 < V < 1.4, Fi- 
nally, for large bias voltage in the range of 1.4 < V < 1.6, 
the width (and area) of the peak gets smaller, resulting in 
decreasing current and the appeareance of negative dif- 
ferential resistance (NDR). Also note that there are weak 
NDR effects in the O.b <V <1.0 region both for the Au2 
and Au3 systems. For H2Pch (OFF-spectrum), the pat- 
tern of IVCs is similar to that of the ON-spectrum, but 
different in magnitude. 



7 



(a) 





A 


K 1.6V 




1.4 V 


A . JK 


11 * 


1.2 V 




A 


l.OV 




0.8 V 








0.4 V 








0.0 V 



(a) 



-2.0 -1.6 -1.2 -0.8 /-0.4 0.0 \Q.4 0.8 1.2 1.6 2.0 



Energy (eVJ\ 



(b) 



S 

1/3 



11 ^^ 


1.6 V 


.A J i 


1.4 V 








l.OV 








1 1 0.6V 


Vr 


0.4 V 

* 


iU 


0.2 V 


— 1 1 1 1 1 1 — jAt 1 — 


0.0 V 



-2.0 -1.6 -1.2 -0.8/ -0.4 0.0 Ch4 0.8 1.2 1.6 2.0 

/ Energy (eV) 



if 



FIG. 5: (Color online) Bias dependent transmission of 
Au(100)-Au3-H2Pc-Au3-Au(100) as a function of energy E. 
The Fermi energy is set to zero; the shaded area indicates the 
bias voltage window, (a) Top panel: ON-state; (b) bottom 
panel: OFF-state. 



The bias voltage drop over Aii(100)-Au3-H2Pc-Au3- 
Au(lOO) in the ON-state in Fig. iTfa) explains much of 
this picture: like for cumulene leads (Fig. |6ja), the main 
bias voltage drop occurs at the left AU3-H2PC contact, 
and most of the remaining bias voltage drop is distributed 
over the bulk of the H2PC molecule, the levels therefore 
basically following the chemical potential of the right lead 
with potential /i = EY-eV/2 decreasing under applied 
bias. 

However, there is in fact a clear difference: close in- 
spection shows that the voltage drop in fact moves, with 
increasing bias across the left benzene ring toward the 
right lead, until around V = 1 it jumps to the right H2PC- 
Aus interface, the chemical potential /i = Ep-{-eV/2 of 
the left lead now extending into the bulk of the H2Pch 
molecule. This is precisely where in the transmission 
spectrum in Fig. [sja) (top) the transmisison peaks be- 
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FIG. 6: (Color online) Electrostatic difference (Hartree) po- 
tential AVh = Vh{V) — Vh{V = 0), showing the bias- voltage 
drop over the junction in the (a) ON and (b) OFF states of 
cumulene-H2Pc-cumulene. For negative bias polarity one gets 
the mirror image of the bias voltage distribution. 



gin following the left lead with upgoing potential. This 
means that the H2Pch molecule now begins to get nega- 
tive charge, pinned by electron transfer from the left lead. 
The picture is also seen in the behaviour of the weak 
LUMO following the charging of the H2Pch molecule un- 
til it crosses the left upgoing leads, gets pinned, and then 
follows the left lead. 

For the H2PCV OFF-case, the transmission spectrum is 
similar to that of the ON-spectrum, but different in mag- 
nitude, the transmission peaks being narrow and less in- 
tense, reflecting much weaker coupling to the leads than 
in the ON-state. A particularly interesting feature is 
the different bias dependence of the leading transmission 
peaks for V > 1, essentially following the chemical po- 
tentials of different leads. This indicates that the HOMO 
and HOMO-1 orbitals at high bias are localized in quite 
different parts of the junction. 

Finally, we conclude that the fundamental reason for 
the ON/OFF IVCs in both cumulene-H2Pch-cumulene 
and Au-H2Pch-Au is again the presence of two interact- 
ing, interfering, conduction channels around the Fermi 
level, being controlled by the tautomer switch. Although 
the bias-dependent charging of the molecular orbitals 
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FIG. 7: (Color online) Electrostatic difference (Hartree) po- 
tential AVh = Vh{V) — Vh{V = 0), showing the bias voltage 
drop over the junction in the (a) ON and (b) OFF states of 
Au(100)-Au3-H2Pc-Au3-Au(100). For negative bias polarity 
one gets the mirror image of the bias voltage distribution. 

controls the positions of the transmission resonances rel- 
ative to the bias voltage window, it is the interference ef- 
fect that physically enhances the difference between the 
ON and OFF state transmissions. This is manifested in 
the Fano-type of line profiles characterizing the transmis- 
sion peaks around the Fermi level. 



IV. COMMENTS ON NDR EFFECTS 

The NDR effect in the case of Au-H2Pch-Au around 
1.6 V bias arises from the decreased intensity (area) of 
transmission peaks at high bias voltage, as shown in Fig. 
[5] The reason is reduced coupling between the Au-leads. 
Since the main conduction peaks follow the chemical po- 
tential of the left lead with upgoing potential, that part 
of the H2PC molecule should be well connected to the left 



lead. It is then reasonable to conclude (and supported by 
transmission eigenchannel analysis) , that the coupling to 
the right lead is reduced. 

Negative differential resistance is of great interest be- 
cause NDR elements can be used for building bi-stable 
devices and highly funct ional components like latches, 
oscillators, memorie^^^ESl ^^^^ ^Qg^^. XOR gate4^^^. As 
a consequence, NDR has been at the focus of a number 
of recent experimental^^ and theoretical^^ papers. 

The NDR effect found by Quek et alW^ is explained in 
terms of gold electrodes supporting molecule-electrode 
contact driven NDR, where the coupling to the rather 
narrow bands of ID Au chains stretched on the surface 
plays an important role. This looks similar to our case 
at high bias. 

The NDR effect investigated by Riberio et alW^ deals 
with quantum transport properties of porphyrin-bridged 
p-n junctions with Si leads, and the NDR effect is con- 
trolled by the doping levels of the leads and the cen- 
tral transition metal atom of the porphyrin. The results 
are explained in terms of bias-induced on-off switching 
of resonant tunneling channels associated with specific 
molecular orbitals. This NDR can then be regarded as 
arising from bias-driven reduced molecule- lead overlap. 
This again looks similar to our case at high bias. 

Finally, the NDR effect studied by Pati et al^ is ex- 
plained in terms of reduced overlap between conduction 
links through the molecule itself. This is in contrast to 
our case where the reduced overlap appears at a contact 
to the molecule, the transmission peaks basically follow- 
ing one lead or the other. In the case of Ref. [60[ the 
peaks do not move with bias, suggesting that the voltage 
drop over the molecule is distributed in a balanced way 
over and around the molecule. 

From an IVC point of view there is no big difference 
between these cases. The differences are primarily visi- 
ble in the T{E^ V) spectrum and reflect the details of the 
bias-dependence of the bias- voltage drop over the junc- 
tions, manifested in strongly non-linear IVCs. With more 
experience it might be possible to relate different types 
of NDR-features to different types of weak-link distri- 
butions and instabilities. It seems likely that voltage- 
or current- driven excitations and instabilities should be 
sensitive to the distribution of strong fields (potential 
drops) over contacts or inside molecules. 

Also note in Fig. [2] the appearance of weak NDR 
around 0.7 V bias in Au(001)-Au2-H2Pc-Au2-Au(001), 
and very weak NDR around 0.9 V bias in Au(OOl)- 
Au3-H2Pc-Au3-Au(001). The latter case can in princi- 
ple be understood from Fig. [S^i: the NDR (or plateau) 
seems to be connected with lineshape variation of the 
HOMO/HOMO-1 transmission peak structure. The 
Fano lineshapes suggest that interference effects again 
may play a role. 
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V. CONCLUDING REMARKS 

We have shown that the two tautomer configurations 
of phtalocyanine (H2PC) connected to conducting leads 
form junctions with widely different conductances, espe- 
cially at low bias. With the H-H pair rotated (hydro- 
gen tautomerization performed) by external means, the 
junction presents a single- molecular switch. This fact, to- 
gether with the experimental demonstration of Liljeroth 
et al^thdit the switch works in principle in STM configu- 
ration (for vertical transport), makes H2PC an interesting 
candidate for a unimolecular electronic switch. 

We find that hydrogen tautomerization affects the elec- 
tronic state of H2PC by switching the order of the molecu- 
lar orbital energies between HOMO and HOMO-1, caus- 
ing a large change in the current. We specifically re- 
veal the influence of molecular charging and molecule- 
electrode coupling on the transport mechanism. In par- 
ticular, we show that the applied bias potential drop can 
shift from one metal-molecule contact to the other at in- 
termediate bias, leading to pronounced features in the 
IVC. We also find that a two-channel interference effect 
is clearly manifested in pronounced Fano shapes charac- 
terizing the transmission peaks close to the Fermi level. 
This leads to the conclusion that the underlying reason 
for the ON/OFF IVCs in the studied configurations is 
the presence of two interacting, interfering, conduction 
channels around the Fermi level, being controlled by the 
tautomer switch. This interference effect is dramatic in 
the sense that it leads to a large ON/OFF ratio over a 
wide bias voltage range. 

Two final comments on the interference effect: (1) 
Since we expect the substrate to play a pronounced role 
stabilizing the two states of the molecule, the question 
arises how this will infiuence the interference, which is 
at the heart of the junction transport. Since the two 
nearly degenerate tautomer states have similar spacial 
localization, we expect the interaction with the substrate 
not to significantly change the level splitting. The levels 
could be shifted relative to the Fermi level, changing the 
conductivity, but that is not really the issue here: the 
ON/OFF switching contrast should remain. (2) The cal- 
culation is done for frozen geometry, taking into account 
redistribution of electronic density on the molecule. The 
shift of the bias voltage drop in Fig. 7 in the range 1- 
1.4V is connected with the steep increase of conductivity 
in Fig. 2. Major electronic density redistribution can be 
expected to lead to conformational changes, and the large 
drop of the conductivity (NDR) in thel.4-1.6V range in 
Fig. 2 certainly needs to be investigated taking into ac- 
count effects of atomic relaxation. But, then again, a 
free- molecule calculation may not be relevant: the atomic 
relaxation will be constrained by the stabilizing effect of 
the environment /substrate. For this reason, we expect 
the calculated NDR effect in Au-H2Pc-Au in Fig. 2 to 
be relevant for the molecule lying fiat on a substrate. 

We also conclude, that the D2/1 molecular symmetry 
of the H2PC with the HOMO and HOMO-1 occupying 



opposite 'arms' with/ without hydrogens is an essential 
component in creating the potential profile across the 
junction. It is important to use this symmetry in order 
to obtain different conductances for different H2PC con- 
figurations. The connection to the electrodes must not 
be to neighbouring arms (in which case the H-H pair ro- 
tation only fiips the mirror symmetry of the junction): 
connecting neigboring arms produces an IVC that is the 
average of the ON- and OFF-states, losing the switch- 
ing effect. Coupling to the electrodes should also respect 
the symmetry by introducing minimal distortion. In our 
study H2PC is physisorbed at the gold electrodes; our at- 
tempts to chemisorb H2PC on gold via two thiol groups 
directly connected to the molecule only resulted in lock- 
ing of the HOMO and HOMO-1 relative positions, thus 
preventing the switching. However, we exemplify with 
cumulene electrodes that it is possible to design linkers 
in the form of molecular ligands to interconnect H2PC to 
external leads. 

We demonstrate, by considering different electrode 
types, that the conduction of the junction is significantly 
affected by the choice of contacts. Both cumulene^'' and 
gold leads are highly conductive, however matching of the 
gold s-states to the H2PC conjugated pi-orbital introduces 
additional backscattering, yielding 4 times less current at 
comparable voltages. In addition, the IVCs of junctions 
with gold leads exhibit pronounced negative differential 
resistance (NDR) at high bias voltage, as well as exam- 
ples of weak NDR at intermediate bias. This, together 
with switching, makes H2PC an interes ting c andidate for 
a multifunctional molecular componentP^'^. 

We realize, that the modelled system with the two 
closely spaced gold contacts is difficult to create in ex- 
periment. Nevertheless we believe that this difliculty is 
an additional good reason for preliminary theoretical as- 
sessment of the setup. Our investigation of H2PC is in- 
spired by the need for switches and non-linear elements 
in molecular electronic circuits. However, the goal of this 
work has been strictly limited to investigating the trans- 
mission properties of the two tautomer states. We have 
not investigated the switching mechanism or stability of 
the OFF and ON states, or the influence of different sub- 
strates. It remains to investigate how the current through 
the junction may switch between OFF and ON states at 
elevated bias, including hysteretic behavior due to relax- 
ation, current-driven fluctuations and thermal noise, as 
well as the influence of substrate surfaces on the conduc- 
tion mechanisms. 
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